Supersymmetry for disordered systems with interaction 
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Considering disordered electron systems we suggest a scheme that allows us to include an electron- 
electron interaction into a supermatrix cr-model. The method is based on replacing the initial model 
of interacting electrons by a fully supersymmetric model. Although this replacement is not exact, it 
is a good approximation for a weak short range interaction and arbitrary disorder. The replacement 
makes the averaging over disorder and further manipulations straightforward and we come to a 
supermatrix cr-model containing an interaction term. The structure of the model is rather similar 
to the replica one, although the interaction term has a different form. We study the model making 
perturbation theory and renormalization group calculations. We check the renormalizability of the 
model in the first loop approximation and in the first order in the interaction. In this limit we 
reproduce the renormalization group equations known from earlier works. We hope that the new 
supermatrix a-model may become a new tool for nonperturbative calculations for disordered systems 
with interaction. 

PACS numbers: 71.23.An, 72.10.Bg 
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I. INTRODUCTION 

Nonlinear a-models serve as an important tool for the 
theoretical description of disordered and chaotic systems. 
These models efficiently describe physics at large dis- 
tances and times, which is the most important hmit in the 
theory of localization, mesoscopic fluctuations, etc. The 
main idea of the approach is to integrate out fast "elec- 
tronic degrees of freedom" and to reduce calculations to 
solving a model containing only slow "diffusion modes" . 
This is the way how nonlinear cr-models are derived. The 
initial model of Ref.^ was not completely satisfactory be- 
cause it contained formally divergent integrals. In order 
to avoid such integrals one should either shift contours 
of integration in a complicated way^ or start with a rep- 
resentation of the Green functions in terms of functional 
integrals over anticommuting (Grassmann) variables^. 

Although these replica cr-models proved to be very con- 
venient for writing perturbation expansions in diffusion 
modes and renormalization group equations, it became 
clear soon that they were not very useful in nonperturba- 
tive calculations. As an attempt to improve the method, 
a supersymmetric version of the field theory was formu- 
lated by Efetov^. The supersymmetry method combines 
both the integration over commuting and anticommuting 
variables and there is no need to use the replica trick. It 
turned out that the supersymmetry method was very effi- 
cient for studying various problems of disorder and chaos 
(for a review, see, e.g., Refj^). 

The next step was to include electron-electron interac- 
tions and this was done by Finkel'stein& who generalized 
the (7-model of Refi^ to interacting systems using the 
replica approach and obtained a set of renormalization 
group equations. It is interesting to mention that, for 
electron systems, the interaction can only be included 



using the representation of Green functions in terms of 
integrals over the anticommuting variables suggested in 
Ref.'^. The alternative representation of Refs.^'^ may 
work for interacting bosons only. Needless to say that 
the CT-niodcl for interacting systems is considerably more 
complicated than that suggested for noninteracting ones. 

One more possibility to derive a cr-model for disordered 
systems is based on the Keldysh technique^. As both the 
replica and supersymmetry approaches, the Keldysh for- 
malism enables one to average over the disorder in the 
beginning of all calculations, which is the main step when 
deriving the cr-models. The special form of the Keldysh 
time contour automatically ensures that the weight de- 
nominator in the representation of the Green functions 
via functional integrals is equal to unity, which is the ba- 
sis for performing the averaging over the disorder. The 
Keldysh formalism for the disorder problems was first 
suggested in Ref.* and later used in Ref;^- for deriva- 
tion of a cr-model for non-interacting systems. More re- 
cently, CT-models based on the Keldysh formalism and the 
fermionic representation of Ref^ were derived for inter- 
acting systemsi2iii. These models, as well as the previous 
ones, are applicable in the diffusive regime. 

Although all the a-models listed above are very simi- 
lar as far as their structure and methods for perturbative 
expansions are concerned, they are quite different with 
respect to nonperturbative calculations. The simplest 
case where this difference can be traced is the problem of 
the level statistics in a disordered metal grain. In the cr- 
model language the level statistics problem is equivalent 
to solving a zero-dimensional cr-model. For the super- 
symmetric cr-model, one can reduce the calculation of 
the level-level correlation function to the computation of 
a definite integral over 2 or 3 variables (depending on 
the symmetries of the model) . The latter can be calcu- 
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lated without big difficulties and one comes to the famous 
Wigner-Dyson statistics (see, e.g., Refii^). The method 
of calculations based on supersymmetry can also be used 
for calculations for random matrices!^. 

The situation is quite different when trying to ap- 
proach these problems with either replica or Keldysh 
CT-models. The first published observation of difhcul- 
ties encountered when calculating a proper integral over 
N X N matrices (N is the replica number must be put to 
zero at the end) was presented in Ref^ii. Computation 
with the Keldysh cr-model is not easier because one has 
time (or energy) as an additional variable and integrals 
that appear for the zero-dimensional tr-model are func- 
tional integrals rather than definite integrals over several 
variables. Recently, several attempts have been under- 
taken in order to improve this situation and to develop 
methods of nonperturbative calculations for noninteract- 
ing replica^^'^^ and Keldysh-^- cr-models. Taking into ac- 
count nontrivial saddle points the authors of these works 
succeeded in reproducing oscillating asymptotics of the 
Wigner-Dyson formulas. The latest achievement in this 
direction is apparently the publicationiSi where the level- 
level correlation function was explicitly derived for the 
unitary ensemble using the replica approach. 

In spite of the rather complicated calculations and 
limited success, the attempts to perform nonperturba- 
tive calculations with the replica and Keldysh cr-models 
were motivated by a probable extension to interacting 
systems. Incorporation of interaction into the supersym- 
metry scheme has not been considered and was generally 
believed to be impossible. This opinion is, of course, 
not groundless. As we have mentioned, in the presence 
of interaction, it becomes crucial whether one considers 
fermions or bosons. The one particle picture that al- 
lows one to consider the fermionic and bosonic represen- 
tations on equal footing is applicable for noninteracting 
particles only. In the supersymmetry method one uses 
integrations over both fermionic and bosonic variables to 
describe, e.g., a system of electrons and it is completely 
unclear how one can apply this method for interacting 
particles. 

There is no canonical way to incorporate interaction 
into a supersymmetric model. In particular, it is hard 
to see how the normalization of the partition function 
to unity could be achieved by adding bosonic variables 
while keeping the full information encoded in the theory. 
In reality, we cannot include an arbitrary interaction into 
the supersymmetric cr-model exactly. Therefore our goal 
is more modest and we restrict our consideration to weak 
interactions. 

Instead of trying to include interaction into the super- 
symmetry scheme we propose an artificial model with 
interaction that is fully supersymmetric and allows the 
standard treatment within the supersymmetry scheme 
(the supersymmetry is violated only by source terms that 
allow us to calculate different correlation functions). We 
justify the usefulness of this model by comparing it with 
the initial electron model and find that, for a weak short 



range interaction and an arbitrary disorder, the models 
are close to each other. The comparison is performed 
by writing diagrams perturbative in the interaction but 
exact in disorder. It is shown that on the level of a 
generalized Hartree-Fock-approximation in the interac- 
tion (but exact in disorder) the supersymmetric model 
we suggest and the initial electron model are equivalent 
to each other. 

Having demonstrated this equivalence we derive the 
cr-model using integration over 4M-component supervec- 
tors {M is the number of Matsubara frequencies, spin 
is included). Subsequently a supersymmetric cr-model 
is derived for systems with unitary symmetry using the 
standard scheme of derivation based on the averaging of 
the disorder, integration over supermatrices Q and cal- 
culating the integrals in the saddle-point approximation, 
which fixes the eigenvalues of Q. In the limit of vanishing 
interaction the resulting model reduces to a form simi- 
lar to the noninteracting cr-model, the main difference 
being the size of the supermatrices that carry two addi- 
tional discrete frequency indices. For noninteracting sys- 
tems the supersymmetry technique proved to be a very 
powerful tool for the analysis of nonperturbative effects 
related to disorder. The interesting question arises how 
a weak interaction influences the physics in this regime. 
The present model is intended to provide a starting point 
for tackling problems of this kind. 

Before addressing new problems, one should make sure 
that well-known perturbative results of the theory of in- 
teracting electrons can be reproduced with the super- 
symmetric cr-model derived. We address this question in 
the present paper by studying the renormalization group 
procedure introduced by Finkel'stein® in the framework 
of the replica cr-model. Of course, as our model is written 
for a weak interaction, we cannot make a comparison for 
an arbitrary strength of the interaction. Therefore, we 
restrict ourselves to writing renormalization group equa- 
tions in the first order in interaction. Under this restric- 
tion we are able to follow the scheme of the renormal- 
ization of Ref^ and demonstrate the renormalizability of 
the model. As a result, we obtain renormalization group 
equations that agree with those written in Refj^ in the 
limit of small effective interaction amplitudes. 

This paper is organized as follows. Introducing the 
partition function of an interacting electron gas with dis- 
order in the Matsubara technique we first consider the 
derivation of a cr-model for the theory without interaction 
(Sec.^. Here we point out some technical differences to 
the standard approach to this problem^. Then we turn to 
the interaction part in Sec. IIIII motivating the introduc- 
tion of the supersymmetric model that is then studied in 
the remainder of the paper. First we present a diagram- 
matic analysis (exact in disorder), by which the model 
may be compared to the exact theory and the range of 
applicability can be assessed. Then, a supersymmetric 
cr-model is derived in Sec. IIVI Contraction rules use- 
ful for perturbation theory are given in Sec. and we 
consider the density-density correlation function as well 
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as Altshuler-Aronov correctionsi^ to conductivity. On 
these examples we show how one can perform perturba- 
tive calculations for the model derived. We demonstrate 
that in certain cases one can include also higher order 
interaction effects. Next we derive the renormalization 
group equations at lowest order in the interaction am- 
plitudes (Sec. IVf|) and also obtain the Altshuler-Aronov 
corrections^^ to the density of states. 



II. SUPERSYMMETRIC MODEL WITHOUT 
INTERACTION 

The conventional supermatrix cr-model for noninter- 
acting particles^ is written for two different frequencies. 
The interaction makes all Matsubara frequencies impor- 
tant and we adjust the supersymmetry formalism to this 
case. In order to understand the scheme better we con- 
sider first the noninteracting case. Of course, the final 
results of calculations with the effective model derived 
here must be the same as those obtained in Reff^ within 
the conventional technique but there are some peculiari- 
ties in the intermediate steps. 

The partition function of a disordered interacting 
fermion gas may be written as a functional integral in 
the imaginary-time formalismSS 



Z = J D{x*,x) e- 



(2.1) 



Here the integration is over Grassmann fields XtT(r, t) and 
X^(r, r). The variables x and x* depend on the position 
r, imaginary time r and the spin index a and obey the 
conditions 

Xair, t) = -Xcrir, T + P), r) = -Xa(i", T + l3). 

(3 = T^^ is the inverse temperature. These constraints 
should be respected when introducing the Fourier trans- 
form and we write 



X(r,n) = VT / dTe'-"-x(r,r), 
Jo 

where w„ = (2n+ l)7rT are fermionic Matsubara frequen- 
cies. 

We decompose the action S into three parts, 

S — Sf -\- Sdis + Sint- 

The first term Sj describes free electrons with a chemical 
potential fi, 

Sf = drdr x*A^.t) {^t + ^- /^^ X<T(r,T), 

where m is the electron mass and summation over the 
spin index is implied. 

The disorder part Sdis describes the coupling to a static 
disorder potential U{t) 

Sdis^ dT(irx*(r,r)f/(r)xa(r,T). 
Jo 



When deriving the cr-model we will average over this po- 
tential assuming that it is Gaussian and (5-correlated 

(C/(r))=0, (C/(r)t/(r'))--^<5(r-r'), (2.2) 

Ztttv 

where v is the bare density of states at the Fermi surface 
and T is the elastic scattering time. 

By 5*0 we denote the non-interacting part of the action 

= 5*/ + Sdis- (2.3) 

The interaction part Sint is chosen in a conventional form 

1 

xT/o(r-r')x/^(r',T)xa(r,T). (2.4) 

where Vo(r — r') is the bare potential of the electron- 
electron interaction and the subscripts a, /3 stand for 
spin components. 

In order to calculate correlation functions one should 
add to the action 5* sources of a proper symmetry. 

The aim of this section is to generalize the supersym- 
metry technique to all Matsubara frequencies and express 
the partition function in terms of a functional integral 
over supervectors. In order to arrive at a supersymmet- 
ric form we should add bosonic variables. Some caution 
is necessary to ensure convergence of the resulting Gaus- 
sian integrals, which is very important for consideration 
of nonperturbative effects. To this end we introduce new 
fermionic variables 



X*n jX^sgn(w„)> 
so that now 5*0 takes the form 



(2.5) 



So = -i'Yl y"^r Xa(r,7i) sgn(-w„) (2.6) 

Xa(r, n). 



-iujn + fi + U{r) 

2m 



To simplify the notation and to arrive at a familiar struc- 
ture we arrange the fermionic fields x(r, n) with different 
Matsubara indices n into vectors x(r) in the following 
way 



X(r) 



x(r,-2) 

X(r,-1) 
X(r,0) 
X(r,l) 
X(r,2) 



V 



(2.7) 



With the help of 



A„fc = Sgn {-UJn)Snk 



(2.8) 
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the conjugate vector is defined as x(r) = x^i")^' 
x(r,n) = x^i"; "•)sgn (— Wn)- We also use the frequency 
matrix Enk = -w„(S„fc. 

With these notations we rewrite the action 5*0 in the 
form 



So 



U(r) 



(2.9) 



where 4 = " ^- . . . 

Of course, when writing integrals over the anticom- 
muting variables x ^ change of variables does not bring 
anything new. However, it is quite important for intro- 
ducing supervectors ip because one gets convergent inte- 
grals over bosonic variables. We write these supervectors 
as follows 



\ 



^(r,-l) 
V'(r,0) 
V'(r,l) 



V 



-0(r,n) 



X(r,7i) 
5(r,n) 



where S (r,n) , S* (r, n) are bosonic variables. Now we 
replace 

X ^ V'j X ^ '4' — V''''^ 

in 5*0 to find the supersymmetric action 



S'o^^i I dri^a(r) iE + ^p + U{r) i^„(r). 



(2.10) 



With this preparation a cr-model can be constructed for 
the noninteracting theory following standard considera- 
tions as described, for example, in Refi^. For simplicity, 
we consider the ensemble with the unitary symmetry, as- 
suming that a small magnetic field suppresses Cooperon 
modes. Then, it is not necessary to further double the 
size of the supervectors. As a consequence, after the 
disorder averaging the resulting quartic term may be de- 
coupled with the help of a 4M x 4M (M is the number 
of Matsubara frequencies) supermatrix fields Q{r) with 
two spin and two Matsubara indices 

^-j^ J dr{^{r)4,{r)Y ^2.11) 
DQ e"^ / dri,{r)Q{r)i,(r)~f^ J dr StrQ^(r)_ 

For the definition of the supertrace operation "Str" 
we refer to Refi^. In contrast to the orthogonal case 
no special caution is necessary for the decoupling, since 
Cooperon modes are frozen. The saddle point equation 
reads 



Qir) 



1 



1 



(V(r)V;(r)) = — 5(r,r), (2.12) 



lE + ^p 



27^W 



5(r,r') = *<5(r-r'),(2.13) 



and is solved by Q = A. In complete analogy to the 
standard case-^ one derives the cr-model 

Z = J DQ e-^ (2.14) 

F = ^ J drStT[D{VQ)^ - iEQ]. (2.15) 

The supermatrix field Q{r) obeys the nonlinearity con- 
straint Q'^{r) = 1 and further Q^r) = KQ{r)K. D = 
VpT/d, where d is the dimension, is the classical diffusion 
coefficient. Here and in the following we adopt the con- 
vention that all internal indices that are not displayed 
explicitly are summed over as part of the supertrace op- 
eration. Repeated indices are also summed over. The 
validity of the model is restricted to the diffusive regime, 
momenta are restricted by the condition ql and tem- 
perature by Tt <^ 1. 

III. INTERACTION PART 

In the previous derivation of the cr-model for noninter- 
acting particles containing all Matsubara frequencies we 
closely followed a corresponding part of the paper^ . Now 
we are to include the interaction. Unfortunately, unlike 
the fermionic replica approach of Ref^, the supersymme- 
try does not allow to include the interaction in a simple 
way and one cannot further follow the scheme of Ref.^. 

Instead of trying to derive a cr-model from the original 
electron model with interaction we follow another root, 
namely, we replace the original electron model with inter- 
action by an artificial supersymmetric one. There are two 
requirements for this new model: it should not be very 
different from the original one and it should be treatable 
within the supersymmetry scheme. We have succeeded 
in constructing such a model that gives the same corre- 
lation functions as the original one in the limit of a weak 
interaction. Due to its supersymmetric form the deriva- 
tion of the cr-model can be performed using the standard 
scheme. 

We start our discussion writing the original interaction 
term Smt , Eq. H2.4|l , in a slightly different form using the 
transformation, Eq. (|2.5|l and the vectors Xj X; Eq. (|2.7|l . 

S^nt = drdv' (x"(r)A^x"(r)) (3.1) 

xyo(r-r') (x^(r')A-V(r')) , 
where the frequency matrix was introduced, 

and a and /3 stand for spin components. 

The integrand in the interaction term, Eq. (|3.1(l . con- 
tains a product of two "scalar products" of the vectors x- 
Using the anticommutation relations for the Grassmann 
variables x we can also rewrite Eq. (I3.1|l in an equivalent 
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form 



xVo{r-r') {x^ (r') A-^x" W) 



(3.2) 



Now the crucial step is to guess an interaction term 
that could be written in terms of the supervectors ip. 
We have understood that for the noninteracting part 5*0, 
Eq. H2.H|I . we could simply substitute the vectors x by 
the supervectors ^ and obtain Sq, Eq. (|2.10|l . It is clear 
that just replacing the vectors x by the supervectors in 
Eq. (|3.1() cannot be the resolution of the problem, since 
in this way one takes into account "Fock" diagrams in 
the first order but cannot obtain the "Hartree" ones. In 
contrast, replacing the vectors x by the supervectors ip in 
Eq. (|3.2(l one can reproduce the "Hartree" type diagrams 
but the "Fock" diagrams are zero. 

Actually, this property is the key to the construction 
of the new supersymmetric model. In order to take into 
account both "Hartree" and "Fock" type diagrams we 
write the new interaction term S'^'^j in the following form 



T 



drdv' 



(3.3) 



[(V-" (r) A^'V-" (r)) Vo (r - r') (r') A'V'' (r') 

- (VJ" (r) A^ij'^ (r')) Vq (r - r') (V^'^ (r') A^^'i/;" (r))] 

The structure of the two terms in the integrand in Eq. 
()3.3(l repeats the structure of Eqs. (|3.1|l and H3.2|l . How- 
ever, due to the supervector structure of -0 these terms 
are not equal to each other, which contrast the case with 
the anticommuting vectors x- Moreover, if ip contained 
only bosonic variables these two terms would cancel each 
other completely. 

Equations (|TTn|) and give the action 



S — Sn -V 



(3.4) 



of the new supersymmetric model we want to use. The 
introduction of this model is the main step toward the 
derivation of the supermatrix cr-model. 

Before starting the derivation of the cr-model it is very 
important to understand which important features of the 
original model are kept in the supersymmetric one and 
what is lost. Of course, the partition function Z' 



Z' = j cxp (- 



S')Dtp 



(3.5) 



of the supersymmetric model is exactly equal to unity 
for any parameters of the Hamiltonian, which is a conse- 
quence of the supersymmetry. Therefore, a comparison 
of the partition functions of the models does not make 
sense. At the same time, the partition function of the 
initial model does not contain any singularities in any 
dimension and this is not a quantity we want to study. 

Interesting properties of the model can be seen in such 
quantities as conductivity, tunnelling density of states, 



etc. These quantities are expressed in terms of one- 
and two-particle Green functions. So, we compare the 
initial electron model with the supersymmetric one by 
writing diagrams describing lowest orders of the pertur- 
bation theory in interaction for Green functions. When 
discussing one particle Green functions we compare the 
self-energies of the Green functions, whereas vertices can 
be calculated for two particle correlation functions. The 
following correlation functions are most interesting for 
us: 

The imaginary time Green function. 



5^^-(r,r',T,T') = 



It can be expressed as 



-(x.(r,T)x:,(r',T'))<, (3.6) 
jD{x*,x) (...)e-^ 



^ D{x\x) e-^ 



(3.7) 



g„„, (r,r',T, t') 
= T^^e,,, (r,r',c^„,c^„) e-^-™-+-"-', (3., 

mil 

and calculated with the supersymmetric model 



Qaa' (r,r', w„i, w„) 



- {xa (r,m) Xcy' (r^^^))5' 



D{^,^) (., 



(3.9) 
r^'. (3.10) 



The correlation function H(r, r', t, t'). 

n(r,r', r,r') = ((p(r, r)p(r', r')))^ (3.11) 
p(r,r) = ^x;(r,T)x.(r,r) (3.12) 

cr 

By ((. . .)) we denote the connected part of the averages, 
{AB) = {{AB)) + {A) {B). We rewrite this correlation 
function in the form 

n(r,r', r,T') = ^ H(r, r', 17^, f],) 6"''''="+''''"', 



kl 



(3.13) 

H(r,r',f]fe,r!,) = - (x(r)A-^x(r) ^(r')AV(r'))s, • 

(3.14) 

The average (. . ) g, was defined in Eq. H3.1()|l above. 

Equation (|3.9|) gives the one particle Green function 
for the supersymmetric model, while Eq. (|3.14l) intro- 
duces the density-density correlation function. In both 
cases the action S' is given by Eqs. H2.10|l and (|3.3|l . 
The presence of the variables x in Eqs. (|3.9() and (|3.14|l 
breaks the supersymmetry and allows us to have nonzero 
results. Expanding the exponentials in and calcu- 
lating the functional integrals over i/i in Eqs. 1)3. 9|l and 
(|3.14l) we obtain contributions that should be compared 
with proper terms of the perturbation theory for the ini- 
tial model. Performing this calculation it is convenient to 
use a diagrammatic language, which will be introduced 
now. 

We have three building blocks: 
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(a) (b) 




(c) (d) 

FIG. 1: The vacuum diagrams at first order in the interaction 
potential. All of them vanish due to supersymmetry. 

(a) The Green function of the nonintcracting theory (ex- 
act in disorder) is denoted by a solid line: 



(b) An interaction line corresponding to the first term 
in Eq. H3.3(l ("Fock interaction") is denoted by a 
single dotted line: 



(c) An interaction line corresponding to the second term 
in Eq. H3.3|l ("Hartree interaction") is denoted by 
a double dashed line: 



Of course, both types of the interaction introduced in 
(b) and (c) correspond to the bare interaction Vq {r — r'). 
The difference between them is related to how they are 
connected to the superfield -0. 

We can check easily doing the perturbation theory that 
the partition function Z' is normalized to unity. This 
property will be used in the next section when perform- 
ing the disorder averaging as for the noninteracting case. 
The normalization is obvious for the noninteracting part. 
On the other hand, for a perturbative expansion in S'j'„j 
the corresponding Gaussian integrals taken with S'q van- 
ish, since the resulting expressions contain products of 
terms of the form Str(G") — 0, where here G denotes 
the average (V'V') taken with S'q. Diagrammatically, in 
the first order in the interaction one finds four diagrams 
depicted in Fig.^ Figures (16) and (Id) can be written 
as a single supertrace, Figs, (la) and (Ic) as the product 
of two supertraces. Closed loops in the diagrams corre- 
spond to supertraces, if we agree that each line of the 
interaction (c) may be used to close a loop. Now we see 
why it is convenient to use the single and double lines 
for the interaction. With this convention all the four 
diagrams in Fig. ^ contain closed "superparticle" lines 
and vanish due to the supersymmetry. According to our 
convention, Figs. (15) and (Id) contain one loop, while 
Figs, (la) and (Ic) contain two loops. 

With this intuition we can discuss interaction correc- 
tions to Green functions depicted in Fig. |21 The first 
[Fig. (2a)] and fourth [Fig. (2c?)] diagrams vanish, while 




(a) (b) (c) (d) 

FIG. 2: Corrections to the Green function at first order in 
the interaction potential. Diagrams (a) and (d) vanish, di- 
agrams (b) and (c) give the Hartree and Fock contribution, 
respectively. 

/••■•o-...^ ,^.o^..^, /-O^^ 



(a) (b) (c) 

FIG. 3: Screening of the interaction lines. These are the 
diagrams at second order in the interaction potential. The 
first and second diagrams are equal to zero, the contribution 
of the third diagram is too large by a factor of 2. 



the second [Fig. (26)] and third [Fig. (2c)] give the Fock 
and Hartree contributions. Note however, that screening 
of the interaction lines inside these diagrams is not cor- 
rectly reproduced. While the first [Fig. (3a)] and second 
[Fig. (36)] diagrams vanish, the third [Fig. (3c)] diagram 
gives twice the desired contribution. The reason is a sim- 
ple combinatorial factor related to the fact that we work 
with two interaction terms. 

This shows that our replacement of the initial electron 
model by the supersymmetric one is not exact and its 
accuracy is restricted by the lowest order in the interac- 
tion for the self-energy. The screening of the interaction 
lines in the self-energy is described by higher orders in 
the interaction and is not given correctly by the super- 
symmetric model. 

At the same time, the density-density correlation func- 
tion n, Eq. H3.14|) . can also be correctly described by the 
supersymmetric model. In addition to the corrections 
to the self-energy of the Green functions the density- 
density correlation function contains vertex corrections 
represented in Fig. 01 In order to write correctly this 
function in terms of the integral over the supervectors 
one should have both supersymmetric and nonsupersym- 
metric vertices, which corresponds to the definition, Eq. 
(|3.14|) . Again, using our two types of interaction lines 
we can check that the ladder diagrams for the density- 
density correlation function are reproduced correctly. Di- 
agrams containing closed supersymmetric loops (the sec- 
ond and the fourth diagrams in the first line and the first 
diagram in the second line) are equal to zero but the re- 
maining diagrams (e.g., the first and the third diagram 
in the first line of Fig. 0] and the second diagram of the 
second line) do the job properly. One can easily extend 
this to higher order ladder diagrams. 

Thus, we see that the replacement of the initial electron 
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FIG. 4: First order diagrams for the density-density correla- 
tion function. The symbol • stands for a source term that 
distinguishes between fermionic and bosonic variables, while 
o symbolizes a supersymmetric source term. Choosing the 
source terms in this way, no artificial contributions are gen- 
erated in this approximation. 



model by the supersymmetric one allows us to reproduce 
reasonably the main order in the interaction. Let us em- 
phasize that the consideration of the perturbation theory 
in the interaction presented in this section was performed 
using Green functions for a given disorder potential. No 
averaging over the disorder or making approximations 
with respect to it was necessary. 

At the same time, our aim is to derive a cr-model, which 
implies an averaging over the disorder. The form of the 
supersymmetric cr-model allows us to perform the aver- 
aging over the disorder in exactly the same way as it has 
been done for the noninteracting case. The scheme of the 
derivation is similar to the one used by Finkel'steinSiSi 
and we can simply follow his procedure. 



IV. 



a-MODEL WITH INTERACTION 



Now we construct a low energy theory for the model 
introduced in the preceding section. After the disorder 
averaging and the decoupling of the quartic term by inte- 
gration over the supermatrix field Q{r) in complete anal- 
ogy with what has been done for the noninteracting case 



J 



one finds 



Z' 



D{iP,iJj)DQ exp ( / drStiQ^ir] 



X exp J dr ■0(r) 
X exp(-S'-„J-!/',V']) 



^E + Cp + —Q{r) 



(4.1) 



where the interaction term S'^„^ is given by Eq. H3.3(l . 

As in Ref A, we should single out in the quartic inter- 
action S'^^f [ip, 4'] pairs slowly varying in space and time. 
However, in contrast to the replica approach, our super- 
symmetric interaction, Eq. H3.^{|l . contains two terms. 
Hence, after singling out the slow pairs we obtain four 
terms instead of the two terms in Ref.^. At the same 
time, one obtains as usual only two interaction ampli- 
tudes Fi and r2 describing the small and large angle 
scattering. These amplitudes can be obtained by the re- 
placement Vb(q) — > Ti/i/ and Vb(pi — P2 ~ 'i) ^2/1^ 
for small momenta q. The quantities Fi and F2 allow to 
make contact with the Fermi liquid theory. In particu- 
lar, the singlet and triplet amplitudes Fg = Fi — and 
Ft = iF2 are related to the Fermi- liquid constants F° 
and F°, 



1 



21 + F^' 



Ft = - 



1 F° 



2 1 + FO ■ 



(We remind, however, that we may work only in the 
limit of small Fi and F2). 

When integrating over the fields -0,-0 we consider the 
simplest possible approximation replacing exp (—Sl^f) by 
exp(— (S'j'nt)) f^or the interaction term, where averaging 
is with respect to the quadratic form in in Eq. (|4.1|l . 
Using the relation {^■ip{r)tp{r')'j = ^(r, r'), where g{r,r') 
is the Green function corresponding to the quadratic in 
-0 form of Eq. (|4.1() , we get 



i^int) — 



^ J (dpi)(dp2)(dq) Str[A^g""(pi - q,pi)] Vo{ci) Str[A-^/^'(p2 + q,P2)] 

I J (dpi)(dp2)(dq) Str[A^g'^"(pi - q, pi)A-^g"^(p2 + q,P2)] Vo{pi - P2 - q) 

^ J (dpi)(dp2)(dq) Str[A^/"(pi - q,pi)] Vo{pi - P2 - q) Str [A-^5"^(p2 + q, P2)] 

^ y'(dpi)(dp2)(dq) Str[A^g""(pi-q,pi)A-^/'5(p2-fq,P2)] Vo{q). 



(4.2) 



From the momentum fiow in the interaction potential one 
reads off processes with small and large angle scatter- 
ing, which should be described by amplitudes Fi and F2, 
respectively. The singling out of the slow modes corre- 



sponds to a restriction of small momenta q and frequen- 
cies. Then, calculating the integral over Q we may use 
the saddle point approximation that gives for the saddle 
point (?(r, r) = 7ri^Q(r). Then, the effective free energy 
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can be written as 

Fint = {S'mt) 

+Str[AJQ""A-^Q'3/3j| 
^{nup^ J ^^^^ {str[A-'Q'5"] Str[A-JQ"^] 

+Str[AJQ'^"A-Jg"'3]} . (4.3) 

Working with the saddle point of the noninteracting the- 
ory is a frequently used approximation in the context 
of cr-models for interacting disordered systema^iiS^. The 
derivation for the remaining parts of the free energy can 
be performed exactly in the same way as in the replica 
approach. As a result, we come to the free energy func- 
tional F [Q] describing the supermatrix cr-model, 

F[Q] ^ ^ J drStT[D{VQy -4EQ] 

+ ^Jdr [T2Q12Q - riQ7iQ] (4.4) 

where 

Qj,Q = 2^r^ (Str[g,^f„JStr[Q^^f„J 

riiafJ 

I 



+Str[Q««„^O^f„J) 5„,+„3,„,+„„ (4.5) 
Q72Q = 2^r^ (Str[Q,^f„JStr[Q^3"„J 

+str[g,^f„^g^^„J) (4.6) 

Again, we note the constraints Q^(r) ~ 1 and Q^(r) = 
KQ{r)K, where K is the matrix 

Matrices K and k act in the advanced/retarded {AR) 
and boson/fermion (BF) sectors, respectively. 

The interaction part of the supermatrix a-model, Eqs. 
(|4.5|) and H4.6|l contains, in contrast to the replica cr- 
model, in total four terms and is somewhat more compli- 
cated. At the same time, its form is quite similar to the 
one in the replica approach. It is important to emphasize 
that the u-model obtained is fully supersymmetric as it 
should be. The supersymmetry can be violated by source 
terms that can be added for generating correlations func- 
tions. 

For actual calculations it can be convenient to express 
the interaction part in a different way, using auxiliary 
fields ipi,4)i. 



F^nt = j dvdv' (^(str[^^(r)Q(r)]Str[^^(r')Q(r')])^^ + (str[^^(r)Q(r)]Str[0^(r')O(r')])^J 

+ jdvdv' (^(str[^^(r)g(r)]Str[^^(r')Q(r')])^^ + (str[0^(r)Q(r)]Str[0^(r')Q(r')])^J. (4.7) 



In this way one can work with expressions including only 
supertraces and not single components. This makes the 
structure more transparent and additionally it has tech- 
nical advantages when working with supermatrices, since 
anticommutativity of the ferniionic components is taken 
care of automatically. The fields Lpt, (pi do not have a non- 
trivial structure in the AR space. Instead, this structure 
is given by 



^(r)=T^(/.,(r,j)A^" 



where j is a frequency index. The same definition is used 
for (pi. The fields fi,f2 have only commuting entries 
and are proportional to the unit matrix, while (jji , (f>2 are 
conventional supermatrices. For convenience, the explicit 
form together with the Gaussian integrals used to define 
the averaging are given in the appendix □ Of particular 



importance are relations giving the averages of the fields 

1(5 

m 

= S"J, ^6{r-r') (4.8) 



^1 (r, (r', -j) + 01 (r, (r', - j) ) 



and 



(r, j)'P2 (r', - j) + 02 (r, j)(t>2 (r', - j) ^ 



aS 



2 6' 



aS ±2_ 

vT 



(4.9) 



where upper indices refer to spin, lower indices to 
fermion-boson space. The importance of these relations 
will become obvious later. 

Using the formalism developed in this section we can 
express the quantities like Green function in coinciding 
points determining the tunnelling density of states or the 
density-density correlation function in terms of a func- 
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tional integral over the supermatrices Q with the free en- 
ergy functional F, Eq. H4.4|l . This can be done using Eq. 
(|4.HI and integrating over ip in Eqs. (|3.9|l and l|3.14(l with 
the effective free energy functional of Eq. H4.ll) contain- 
ing both Q and ip. The following simple formulas are ob- 
tained when averaging the interaction part and the pre- 
exponential factors in Eqs. (|3.9|l and (|3.14|) separately. 
Furthermore the saddle point equation Q{r) = ni'g{r,r) 
(j??T^ is used, 

g,,, {v,v,u;m,iun) ^ '-^ J DQ Str[fc+Q^';(r)] e"^. 

(4.10) 

For the density-density correlation we consider only one 
out of the two possible pairings for the pre-exponential 
factors in Eq. (|3.9(l . since we are interested in points r, r' 
that are far apart and g{r, r') falls off at distances |r — r'| 
of the order of the mean free path, 



= -{-Kvf j ZJQ Str[fc+A-'=Q(r)]Str[A'Q(r') 



(4.11) 



Equations 14. 411 . (|4.5I) . H4.6|l . (|4.10l) . and H4.11|l are the 

final results of the derivation of the a-model with inter- 
action. Any physical quantity we are interested in can be 
computed, at least in principle, using these formulas. Of 
course, our ultimate goal is to compute some new quan- 
tities in the nonperturbative regime. However, in order 
to see how the model works we should understand first 
how it can help to reproduce known perturbative results. 
Such calculations will be carried out in the next sections. 



V. PERTURBATION THEORY IN DIFFUSION 
MODES 

Disordered systems with interaction can be studied 
using a conventional diagrammatic technique (see, e.g., 
Rcf.^^). In low dimensions, the most important con- 
tributions come from a certain class of diagrams called 
"cooperons" and "diffusons" . Although these diagrams 
are important even without any interaction, the latter 
leads to new contributions. The cr- model formalism al- 
lows one to integrate over "electron degrees of freedom" 
in the beginning of all calculations and reduce them to 
study of the low lying diffusion modes. In the language of 
the (T-model summation of the diffusons and cooperons is 
equivalent to an expansion in small fluctuations near the 
ground state of the free energy functional F, Eq. H4.4|l . 
The minimum of F is achieved setting Q = A, where A 
is defined in Eq. H2.8|l . 

For a perturbative expansion the matrix Q may be 
parametrized in the vicinity of A as 



where P = KP'^K 
P = 
B = 



B 

kB"! 



(5.2) 



a la 
t ih 



P 



kB'^ 



Following the structure of "0 all matrices in Ai?-space are 
now arranged in block form as 

M(fc,Z)| 

The perturbation theory in the diffusion modes is car- 
ried out by making expansion in the supermatrices P. 
Making expansion in P we represent first the free energy 
functional F in a form of a series in P. As the zero ap- 
proximation we take the noninteracting quadratic part 

Fo=Tny J dr Str[i:>VF(r)VP(r) + 2EAP^{r)] . (5.4) 

The rest of the functional consists of quadratic terms 
coming from the interaction, Eas. (|4.5|) and H4.6|) . and 
higher order in P terms coming from both noninteracting 
and interacting parts of the functional F. For explicit cal- 
culations we must calculate averages of different powers 
of P with the bare functional Fg, Eq. (|5.4|) . Calculating 
Gaussian integrals it is not difficult to derive the follow- 
ing contraction rules for arbitrary supermatrices M, N 
(we remind the reader that the unitary case is consid- 
ered) 

(Str [Mk'i'Pikir)] Str [N,,,, Pj,{r')] )^ 
= Str[M^,N^^^,] -l-I?(r-r',fc-0 

xSjkSue{l,l')e{k,k') (5.5) 

and 

(Str [AU^,P^k{r)Nk'i'PH{r')])^ 
StrMey StrNk'i' -StrAUyA StrNk'i'A 

X (r - r',j - k) 6kiS,Mk, k')e{j,j'). (5.6) 

The subscript "0" will be omitted in the following. We 
denote = (M- AAfA)/2, and Af H = {M + AMA)/2. 
The diffusion propagator T) is given by 



V{Y~r',l-k) = j (dq) 



piq(r-r'; 



^q2 + |r!i-fc| 



(5.7) 



and the symbol Q{k,l) stands for 
e{k,l) - 



A{i + iP){i ~ ipy 



(5.1) 



1 for A: > 0, Z > or fc < 0, ^ < 
otherwise 



Correlation functions can be calculated expanding 
Eqs. H4.10|l and 14.11|) ) in P and integrating over this 
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variable with the help of the contraction rules, Eqs. (|5.5() 
and H5.6|l . A few elementary diagrams for the density- 
density correlation function are depicted in Fig. 0] As 
discussed in Sec. IIIII one can read off from the diagrams 
in Fig. 2]that only one vertex term should distinguish be- 
tween the fermionic and the bosonic sector, i.e., contain 
only fermionic variables, while the other remains super- 
symmetric. In this way no artificial contributions arise 
from the diagrams shown in Fig. 2| The same argument 
can be generalized straightforwardly to the ladder-type 
diagrams considered below. Using two vertex terms of 
the same kind would not lead to correct results. In this 
case either all diagrams vanish identically if both vertex 
terms are supersymmetric, or all diagrams in Fig. 01 give 
a contribution, if both vertex terms contain fermionic 
variables only. 

Calculating the density-density correlation function 
n(q, i7) we are interested in small frequencies and mo- 
menta. In the static limit it is related to the compress- 
ibility d^n according to 

dnTi — lim n(q, 0). 

For systems without interaction the relation d^n = 2v 
holds. One can calculate 11 (q, fife) from the supersym- 
metric model representing this function as 



n(q,17ft 



dn 



1 



rn(q,f^O, 



(5.8) 



d^i (i + Ff)2 

(F^F^)~r^ / ^(-q,-fc)n(q,r!fc)i?(q,fc) .(5.9) 

7, J a 



k 

The terms and are given by 

Fi) = ini' / dr StrliJi 



(5.10) 



F^ = mu j drStr[k+ipQ], (5.11) 



where fc+ = (1 + k)/2. The factor 1/(1 + F°)2 is due to 
vertex corrections taken from Fermi-liquid theory. It is 
related to the amplitudes F^ as 1/(1 + F^) = 1 - 2ri 
and also determines d^n = 2i// (1 -I- F°). We consider the 
Fermi-liquid approximationSi, for which the calculation 
amounts to the ladder summation depicted in Fig. |S1 In 
the zero order one should expand Q '--^ A -f- 2iAP in the 
interaction terms, Eqs. (|4.5|1 and 1)4. 6|) . in the free en- 
ergy F and in the pre-exponential, Eq. (|5.9|) . In both 
the cases the saddle point A does not contribute. Only 
the singlet interaction amplitude F^ enters and one must 
sum a geometric series that arises after expanding in the 
interaction terms written in the quadratic approximation 
in P. As a result, we obtain 



n(q,l]fe) = -2iy 



i?q2 + z(°)|r!fc| 



where zl^"^ = 1-1- 2Fi — Using this relation and Eq. 





(DO) 




FIG. 5: Ladder summation for the density-density correlation 
function as discussed in the text. 



one comes to the result 



n(q,f^fc) = 9^7i_-^ 



where DpL ~ {^^Ff)D. The diffusion coefficient, which 
enters the cr-model, can be related to conductivity with 
the help of the equation of continuity 



CT = lim lim n(p, 51) 
This leads to the Einstein relation 



(5.12) 



a = e^duuDpL- 



to hold also after 
below) and in this 
conductivity can be 



This relation can be seen 
renormalizatioK^ (see Sec. I VII 
way interaction corrections to 
obtained. 

In the remaining part of this section, in contrast, we 
describe simple perturbation theory for finding the first 
order interaction correction to conductivity from the a- 
model. A general expression for the calculation of the 
linear response conductivity can be given with the help 
of the original fermionic partition function Z in Eq. (|2.1|1 . 



K{v, r, r',r') 



/^(q^o,ao 



(5.13) 



Z[A] 



M^(r,r)M^(r',T') Z[0] 
I if(r,r,0,0) e-^i"" e*""^ (5.15) 

J Y.T 



A=0 

(5.14) 



Analytic continuation is performed from positive Mat- 
subara frequencies. No summation over /i is implied. 
One comes to Z[A\ after the substitution p ^ p — A 
in Z = Z[Q\. The electron charge e appears in the pref- 
actor. For the cr-model the corresponding replacement is 
VQ — > VQ— ?[-A, Q] leading to a change in the free energy 
F ~ F[0\ ^ F[A\. A \s chosen supersymmetric because 
we should treat fermions and bosons on equal footing for 
the derivation of the cr-model. In order to obtain the ana- 
log of K, however, it is necessary to break supersymme- 
try. To this end one should expand Z'[A\ = J DQ e~^["^l 
up to second order in A and multiply one of the fields A 
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by the projector onto the fermionic sector k+. Since we 
are interested in the leading order correction in the inter- 
action we also expand up to first order in Fint- In order 
to express the result of this procedure in a compact form 
we define 



F 



(1)+ 



A 
(2) + 



■KVlD 

2 

■KviD 



drStr[VQ(r)[A(r),Q(r)]], (5.16) 

dr Str[vg(r)[fc+A(r),Q(r)]](5.17) 

drStr[[fc+A(r),Q(r)][A(r),Q(r)]]. 

(5.18) 



We can now write the analog of K for the supersymmetric 
model, 



- F 



6A,,{r, r)M,,(r',T') 
1 



p(2)+p 



Fa Fi, 



A=0 



4 

E 

i=l 



(5.19) 



Averaging is with respect to the noninteracting part of 
the free energy. We expand the matrix Q in terms of P 
and use Fq (|5.4|I for averaging. Classical conductivity can 
be found at the level of the saddle point Q = A. The only 
contribution comes from K'^ and one finds ctq ~ 2e^vD, 
where the factor of 2 is due to the spin degree of freedom. 
Next we turn to interaction corrections. At order K'^ 
vanishes in the limit q 0. There are two contributions 

K'^{vy,T,T') (5.20) 

~ 2{^ufDdA,dA'^j^ (Str[fc+[Ar^,APrJ[4r,,APrJ] 

X ^(-l)^-(str[4.^/Pr.]Str[^^.^^APr3]) ) 
+ 2{^vfDdA,dA'J^ (Str[A;+[Ar^,A][Ar^,Ap2]] 



X ^ (~l)^-(str[4^/Pr,]Str[^^.^^APr3] 
i=i,2 



After further evaluation one finds a cancellation and thus 
there is no correction including only two diffusons. Con- 
sidering now corrections with three diffusion modes the 
only nonvanishing contributions in the limit q — > can 
be seen to originate from K'^ , where all Q matrices should 
be replaced by 2iAP. At order P^ we therefore have 



Kf\v,T,r\T')^Ain,y)^D^dA,A 



(5.21) 



^ StrAVPrJA:+Ar^,APrJ StrAVPrJ^r.,APrJ 



^ (-1)^- (Str[0^.^^ AP,,] Str[^^. ^^AP^J ) ) 

J = l,2 / 

= 2(. . .)(V1V2)(13)(24) + 4(. . .)(V12)(13)(V24) (5.22) 
+ 2(. . ■)(V13)(12)(V24) 

In the last two lines (|5.22|) the relevant contractions and 
their multiplicities have been written. (V13) symbolizes 
a contraction of VP(ri) and P(r3), etc. After taking the 
limit q — ^ one finds 



Ki{(i^o,n„ > 0) 
lep* 



(5.23) 



T 



d 



(Fi - 2r2) J (dp) Dp^ X 



m>n 



Analytic continuation for the conductivity gives 



6a{uj -> 0) 



2iao 
irdv 



(Fi - 2T2) 



00 Q 

oz 



(dq) 



Dq 



—00 
2 



z coth 



{D(Y — iz) 



Pz 
~2 

(5.24) 



In particular, for two dimensions and after replacing 
r, Vi one comes to the well known resultiS 



(5o- = cr - CTo = eV(27r2)(yi - 2V2)ln(Tr) 



(5.25) 



This agreement demonstrates that the supersymmetric 
cr-model we have derived reproduces correctly the known 
quantum corrections for a weak interaction. However, 
this calculation cannot be considered as the crucial check 
of the consistency of the cr-model because the first order 
correction is not very sensitive to the details of the struc- 
ture of the cr-model. Therefore, in the next section we 
proceed further demonstrating the renormalizability of 
the model and deriving renormalization group equations. 



VI. RENORMALIZATION GROUP 

The renormalization group (RG) procedure in the form 
introduced by Finkel'steinSi will now be applied to the 
supersymmetric model. Finkel'stein considered the prob- 
lem of an interacting electron gas in the diffusive regime. 
Based on the replica approach he developed an appro- 
priate cr-model description and demonstrated the renor- 
malizability up to the first order in t = Xj (A.-k^vD') and 
in all orders in the interaction amplitudes. The scaling 
behavior was studied while lowering the temperature of 
the system. It turned out to be possible to make contact 
to the Fermi-liquid theory. 

Typically, the (Fermi-liquid) interaction amplitudes 
are not small and in the process of renormalization they 
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can even diverg o^i^^ . In the original work^ the case of 
long range Coulomb interaction was considered and the 
Cooperon modes were suppressed (e.g., by a small mag- 
netic field). Appropriate scaling equations for this case 
were derived. Later on the approach has been general- 
ized to many different situations (for a review see Refsi^ 
andSS). By a comparison to the field theoretic solu- 
tion the scaling equations can also be written starting 
from conventional perturbation theory^- and for some 
calculations the diagrammatics can also be conveniently 
uscd^^(for a review see, e.g., Refi^). The authors of 
Ref.^- also considered the case of short range interac- 
tions. A recent discussion in the context of a Keldysh 
CT-niodcl can be found in Ref^. 

The renormalization scheme for the interacting dis- 
ordered electron gas is more complicated than for the 
noninteracting case. "Infrared" singularities arise in the 
theory due to both the frequency and momentum inte- 
gration. It turns out that in addition to the diffusion 
coefficient and the interaction amplitudes a "charge" z 
multiplying the energy matrix must be introduced. It de- 
scribes the relative change of momentum and frequency 
scales. Our aim is to show that the supersymmetric 
model introduced above is capable of reproducing the 
correct RG equations in the limit of weak short range in- 
teractions, i.e., in the first order in t and to leading order 
in the interaction amplitudes. Therefore the procedure 
is used as a test for the new model and it is the scheme 
introduced by Finkel'stein that serves this purpose best. 
As a by-product, the Altshuler-Aronov corrections^^ to 
conductivity and the density of states can also be ob- 
tained in this way. 



A. Renormalization group procedure 

In this section we work in two spatial dimensions and 
rewrite the model of Eq. H4.4|) in the following form: 

F\Q\ = ^ J dr Str[i?(VQ(r))2 - 4zEQ{r)] (6.1) 



2 



dvdv' ^Str[^(r)Q(r)] Str[02(r')Q(r')] 
2 , drdr'(str[^(r)Q(r)] Str[0i(r')Q(r')] 



In order to reduce the number of terms we use the 
fields (pi — (pi -\- ^pi together with the symbolic notation 

(^(i>i4>i^ — {(pifki) + iVi^i) ■ The "charge" z has been intro- 
duced as discussed before. We start the renormalization 
group procedure separating slow and fast modes 

Q^UQoU, Qo^UqAUq, UU = UqUo^1. 



Fluctuations described by the supermatrix Uq are fast 
while those described by U are slow. As a result of this 
separation we write the free energy F [Q] in the form 



i^[Q]-F(o)+i^,„t, where 



= -J- I dr Sti[D{VQof + 2D<S>[Qo,VQo] 



+D[Qo,^^ -AzEUQaU], (6.2) 

j drdr' ((Str[0irC/rQorf/r] 



X StrUi, 



,U,,Qor'UA] )^ - (01 ^ 02) j ,(6.3) 



and * = UVU. 

Next, we parametrize fast modes Qq in the following 
way: 

Qo = A{1 + iP){l - iP)-\ {P,A} = 0. 

One must specify the precise momentum range related to 
the fast and slow modes and additionally account for the 
frequency dependence. This can be done in the following 
way: 

1. Frequencies in the interval Ar^^ < |(jj„| < t^^, < 
A < 1 and momenta in the shell At^^ < Dk^ / z < 

are referred to as fast. 

2. If slow variables U have at least one Matsubara 
index corresponding to the fast frequencies, they 
must be set equal to 1. 

3. If fast variables P do not have at least one fast 
frequency or fast momentum, they vanish. 

The separation of modes in the energy part requires 
some care. The energy matrix must be split into a fast 
and a slow part: E = + E'^. Then 



drStT[EUAP^U] = J drStrlE^UAP^U] 

+ J dr Str[£;-^ AP2]. (6.4) 



Using the parametrization introduced above the free en- 
ergy splits into several parts. First, the original free en- 
ergy is recovered with Q replaced by Q = UAU and with 
the slow energy matrix i?*. In Finkel'steins scheme it is 
sufficient to keep terms of the second order in P for the 
remaining parts. As a result one finds 



(6.5) 



Fi 
F2 
Fe 



TTv j drStr[D\IP\7P + 2zEfAP'^], (6.6) 
2ttvD j drStr[^P,VP]], (6.7) 
-2t:vD j dr Str[($A)2p2 _^ (AP$)2], (6.8) 
2'Kvz j drSiv[E''UAP'^U]. (6.9) 
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The interaction part Fint = Fint,i +^int.2 can be written 
as follows: 

=2(7rzy)2y" drdr' ((Sir\^U,kPrU,] 
xStr[^C/r.APr'C/r'])^ - (^1 ^ h)^ , (6.10) 
F,nt.2 = nT^vf j drdr' ((str[^Qr] 

xStr[0ir,C/r.AP3C7r'])^^ - ((^1 ^ h)^ . (6.11) 

The fast modes can be integrated out in the Gaussian 
approximation and an effective free energy F is generated 
in this way 



F[Q] = - In DPe 



~Fi—F2 — FE — Fint q~Fq 



'F[Q]. 
(6.12) 



The functional F is the appropriate free energy for the 
slow modes, taking into account the influence of the fast 
modes. It turns out that this influence can effectively 
be expressed by a change of D,z,Ti, and r2 with the 
scale determined by A and this is what demonstrates the 
renormalizability. From now on it will be clear from the 
context whether we refer to slow Q or fast Q, so that we 
can simply use a common symbol Q. 



B. Diffusion coefficient 



Here we calculate the corrections to the diffusive part 



Fr. 



Str [D{VQf 



The energy part does not contribute and the renormal- 
ization of the diffusive part can be written in the form 



5Fd = (F„ 



1 



{{FlFnt)) - {{F2F.,nt)) + 2 {{F^Frr^t)) ■ 

(6.13) 



By ((. . .)) we denote the connected part of the averages. 

The term Fint itself does not contain any gradients. It 
is therefore not obvious, how {Fmt) can contribute to the 
renormalization of the diffusion coefficient. Therefore we 
discuss this calculation in some detail. 

According to the contraction rule 1)5. 6|) we have 
(PP) =0. As a direct consequence of this property the 
average (Pmt,2) vanishes. For the part Pi„t,i we should 
distinguish between several possible cases. 

(a) The matrices P can have two fast frequency in- 
dices. Then the t/-modes completely vanish from 
the expression. It is quadratic in P and therefore 
contributes together with Fq to the quadratic form 
in P. This leads to the dressing of the diffusion 



modes in Finkel'stein's approach. In our model 
these higher order interaction effects can only be 
taken into account in special cases (as in Sec. IvTl. 
while in general we have to restrict ourselves to the 
lowest order in the interaction amplitudes for the 
renormalization. 

(b) The matrices P have two slow frequency indices, 
but the momentum is fast. It turns out that in- 
tegrating over P the term is quadratic in the slow 
fields Q and contributes to the renormalization of 
the interaction part (see Sec. IVI l51 below). 

(c) The matrices P have one fast and one slow fre- 
quency index. The diffusion propagator, arising af- 
ter calculating Gaussian integrals in P can be ex- 
panded in slow momenta, which leads to a correc- 
tion to the diffusion coefficient. If it is expanded in 
small frequencies, it contributes to the renormaliza- 
tion of z. The latter case becomes important only 
if the diffusion modes are dressed and we refer to 
RefiSi for details. The former case will be discussed 
next. 



We consider the case (c). 
modes one obtains 



After integration over P- 



{F,nt) = -2nvT^ J drdr' V{r~r',k-l) 



xStr 



4>2 



X (A^^[/r),,(t/r'A-^^)^, 



,(6.14) 



where fc is a fast frequency index while I is slow. In 
Sec. lIVl it has been stressed that Eqs. (|4.8|l and H4.9|l are 
crucial for the model. It is in fact very important for the 

consistency of the theory, that ^(^0^ cc 1. Integrating 
over (p we find 

(P„,t) = -27rT (Fi - 2r2) J (dp)(dq) I5(p + q,k-l) 



xStr 



q,m/ ^ — q,/m 



+ t^q,m;t^-q,im0-fe-l,. 



(6.15) 



Here p is a fast momentum while q is slow. Keeping 
only fast k and p in the diffusion propagator the U de- 
pendence would drop completely. Therefore one expands 
the diffusion propagator up to second order in q. Using 
the identity Str[VC/VC/] = — Str[$^] and approximating 
the frequency summation by an integral, we come to the 
following expression 

(fmt) = (Fi - 2r2) J dr Str [<P^] J (dp) 

X J dn{2D'^p'^V^(j),n) - DV'^{p,Q)) . (6.16) 
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The momenta and frequencies in the region of the inte- 
gration in Eq. (|6.16|) are large and correspond to the 
fast variables. We use the symbol /' to indicate that 
integration is carried out over positive frequencies only. 

The remaining terms in Eq. (|6.13|) can be considered 
in an analogous way. The results are 

- {{FiF.r^t)) = 2(ri - 2T2) J dr Str[$ll2] J (dp) 

X J dn{DV'^{p,n) ^2D^p^V^{p,Q)) , (6.17) 

- {{F2F,„t)) = -(Fi - 2T2) J dr Str ^ 

X J (dp) J dn Dv^{p,n), 

2T2) J dr Str [$"2 
X J{dp)Jdn D2p2 v^{p,n). 



$l|2 _ $-L2 

(6.18) 



(6.19) 



Following Finkel'stein, terms that are parametrically 
small in tDq^ / (zAr^^) <C 1, where q is a typical slow mo- 
mentum, have not been taken into account. Such terms 
can arise from expressions in which all frequency indices 
are fixed to be slow by <I>-modes2i. A comparison of the 
combined contribution of Eqs. (|6.16fl^T^ to the original 
diffusive part Fjj gives 



SD = -- 



1 



-Ji 



(27r)2zy 

J^-Mr.-2T2)Jidp)J'dn ^^^,^^^^3 



(6.20) 
(6.21) 



C. Effective charge z 

As we have mentioned, in the cr-model for interacting 
systems a new effective charge z appears as a coefficient 
in the original tr-model. 



Str[£;Q]. 



Expanding this term in P we write the relevant term of 
the second order in P as 

Fe = 2ttvz 

It is understood that the energy matrix E contains only 
slow indices here. In our case the correction to z origi- 
nates only from one term 



5F, 



{{FeF,,u)) 



Working with dressed interaction lines as in Ref.^^, one 
more contribution from (Fint) should be taken into ac- 
count. The final result, however, is the same in both 



cases and we find 

- {{FEF,„t)) ~ ziVi - 2r2) 



drStT[EQ] J{dp)J dnv'^{p,n). (6.22) 



Equations H4.8(l and H4.9(l were again crucial to arrive at 
this result. The relation dn'D{p,^}) = —z'D'^{p,^l) can 
be used now to obtain 

Sz^ — f (dp) (El - 2r2) [Dip, f^)][^:^;-i . (6.23) 

Equation (|6.23() gives the logarithmic correction to the 
effective charge that will be used for writing renormal- 
ization group equations. 



D. Interaction part 

The leading order corrections in Fi, F2 contributing 
to the renormalization of the interaction part come from 
the average (Fint)- We will see that it is very important 
here to work with both fields and ip. As pointed out 
above in Sec. IVI Bl for this term both frequencies of P 
should be slow, while the momentum is to be fast. After 
integration over P we obtain 



{F^nt)^Y /(^P)^(P'O) 



X drdr' (Str[Qr'^lrQr '^lr' - Qr'p2rQr 4l2r'] ) . . 

- j {dp) V{p,Q) j dr (Fi Q72Q - F2Q71Q) • (6.24) 

We write I'(p, 0) here to symbolize that only the momen- 
tum should be integrated, since frequencies of P-modes 
are slow. As a result, the corrections to Fi and F2 are 



<5ri 



1 



(dp)I?(p,0), 
5T2 = — Fi fidp)V{p,0). 

TTV J 



(6.25) 
(6.26) 



This agrees with the corresponding result obtained from 
the replica model^i, since the dressing of the diffusion 
mode is not effective here. We would like to stress that 
the correction to Fi is proportional to F2 and vice versa. 
From the technical point of view it is of particular im- 
portance here that the terms QjiQ and Q^2Q in the 
supersymmetric model of Eq. I|4.4|) include both the sin- 
gle supertrace and the product of two supertraces. Only 
in this way the correct structure can be reproduced in 
the calculation of (Fint)- 



E. RG equations 

From Eqs. (jOn|) . and (j^:^ appropriate 

renormalization group equations may be derived. They 
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take the form 



1 dt 



1 dz 

1 dVi 
IdlnX-^ 
1 dr2 
t dlnA-i 



-(ri-2r2), 

z 

-(ri-2r2), 



(6.27) 



These equations are valid in the first order in t = 
{{2Tr)'^i'D)~-^ and in the leading order in the interaction 
amplitudes for Fj <C z. In the first order in ri^2 they 
coincide with the corresponding limit of the renormaliza- 
tion group equations obtained in Refsi^ andii for short 
range interactions, where terms of all orders in the inter- 
action amplitudes were included. As our supersymmetric 
model, Eqs. (|2.1Q(I and (j^J, can serve as a good approx- 
imation for the initial electron model only for a weak in- 
teraction, taking into account the terms of higher orders 
in interaction would lead to an overestimation. There- 
fore, we conclude that in the region of applicability the 
supersymmetric cr-model gives the correct renormaliza- 
tion group equations. 

From these equations one can see that the combination 
zi = z — 2Ti + T2 is invariant under renormalization. 



dzi 



d\nX- 



rflnA 



— (Z-2F1+F2) =0 



(6.28) 



This invariance is not incidental but protected by a Ward 
identity and therefore its validity is not restricted to the 
approximations studied here. For more details we refer to 
RefsiS'SiSi, andii. Its importance for the consistency of 
the entire scheme will become obvious in the next section, 
when we discuss once more the density-density correla- 
tion function. 



F. Conductivity 

During the process of the renormalization the mean- 
ing of the coefficient D appearing in the free energy and 
its relation to conductivity remains not very well under- 
stood. In order to clarify its meaning one can study 
the density-density correlation functionS*2iiSi2^. In or- 
der to calculate the density-density correlation functions 
one can add to the free energy functional a source term 
of the form = inv J dr StT[dQ], Eq. ifTTUIl . Then 
corrections to this term can be calculated. In our ap- 
proximation one should consider the average 

This expression covers vertex corrections as well as cor- 
rections to the single particle Green function. These con- 
tributions exactly cancel, so that F^, Eq. H5.10|l is un- 
changed. [This is also true in all orders in F/2i]. We 
would like to emphasize that one may not use F^, Eq. 



(|5.11(l . in this procedure because this would lead to spu- 
rious contributions. 

Now the dynamic part of the polarization can 
be calculated from the renormalized cr-model as in 
Sec. It is generally expected that the thermody- 
namic density of states does not develop logarithmic 
singularitiesSi2ii22422i. Then we come to the result 



n(q,r!fc) 



-2i/ 



Dce + zi\nk\ 



where D is the renormalized diffusion coefhcient. The 
coefficient zi is a constant under renormalization as we 
have seen above, Eq. H6.28|l . and using Eq. H5.8() we find 
again 



n(q, f^fe) = d^n 



De = (1 + Fq)D. With the help of the equation of con- 
tinuity (|5.12|l one recovers the Einstein relation 

(T — e^df^nDg. 

So, this relation remains valid and one can determine the 
conductivity from the renormalized diffusion coefficient. 
As an example, one can make contact to the Altshuler- 
Aronov corrections to the conductivity for the case of 
short range interactions. To this end one should replace 
z by its bare value z = 1 as well as F^ — > 1^ in Eq. 16.27|l . 
where Vi = i^Vo{0) and V2 = i^Vo(p — p') (the bar means 
averaging over the Fermi surface) . As a result one comes 
to the well-known formula, Eq. H5.25|! . 



G. Density of states 

The main purpose of this section is to obtain an expres- 
sion for interaction corrections to the density of states 
from the supersymmetric cr-model derived. The single 
particle density of states (per spin degree of freedom) 
can be written as 



Ke) 



j- y'(dp)3G« (p,e) 



No summation over spin index cr is implied. In order to 
calculate the density of states in the Matsubara formal- 
ism it is convenient to first define the function 



j>(w„) 



-]^Jdr GaAr, r, r,0) e^""^ 



where G denotes the Matsubara Green function of the 
system. It is sufficient to know i'(w„) for positive fre- 
quencies. The density of states v can be obtained from 
I' by an analytic continuation. 

In our formalism we obtain (compare Eq. I3.9|l 

DM = -Zi. {Str[k+QZn{r)]) , (6.29) 
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assuming that the saddle point approximation has been 
employed. Again, we separate fast and slow modes in Q 
and expand Qo up to second order in P. We also expand 
exp(— i^) up to first order in Pint- We should consider 
the correction 



^;> = zi.((Str[/c+Q,7„(r)]F„,,)) 



(6.30) 



The calculation is very similar to that of z and one ob- 
tains the correction to the density of states in the form 



t 



5u = -y-{Ti - 2r2)lnA-i. 



(6.31) 



Equation H6.31|l allows to write a renormalization 
group equation for the density of states 



1 dh 



V dhiX 



-(Fi - 2T2). 



(6.32) 



The validity is, of course, restricted to small effective 
interaction amplitudes Ti/z <C 1. As in Sec. IVI Fl one 
can recover the Altshuler-Aronov correction for the short 
range case that can be written in the standard form 



5v/v ^ -t{Vi - 2V2) In(TT). 



(6.33) 



The same formulas Eqs. H6.31|l - H6.33|) for the density of 
states would be obtained with Replica or Keldysh models 
in the limit of weak interactions. 



H. Discussion of the renormalization group 

The renormalization analysis carried out in this section 
allows us to come to the conclusion that the supersym- 
metric cr-model is renormalizable at least in the first order 
in the interaction amplitudes and inverse conductance 
t — {{2tt)'^i'D)^^ . The structure of the interaction term, 
Eqs. (|4.5I) and H4.6|l is crucial for the renormalizability 
and therefore the renormalization group treatment can 
serve as a very good check of the cr-model. The renor- 
malization group equations derived above (Eq. (|6.27|l ^ 
agree with those given in Refsiii*2iiSi, and^^ in the lead- 
ing order in Ti and for <C 0. 

We do not present here the solution of these equations 
referring the reader to KeisM'^i^, and^^ where solutions 
of more general equations have been discussed for arbi- 
trary values of the amplitudes F^. Unfortunately, even 
these more general equations do not allow for definite 
conclusions in two dimensions, since under renormaliza- 
tion the quotient F2/2; diverges. This leads to a diver- 
gence of the spin susceptibility, which has been inter- 
preted as some sort of ferromagnetic instabilitjiSi. How- 
ever, it should be stressed once again that the renormal- 
ization group calculations were performed here only as a 
check of the new cr-model. For the more complete study 
of the scaling behavior in 2D replica or Keldysh models 
are more convenient. The remarkable point is that the 
nontrivial set of Eqs. H6.27|l has been obtained from a 
model based entirely on the supersymmetry method. 



VII. DISCUSSION 

In this paper we constructed a supersymmetric cr- 
model for disordered fermion systems with interaction. 
Instead of trying to combine the supersymmetry tech- 
nique with an electron-electron interaction for the initial 
electron model we proposed an artificial fully supersym- 
metric model with interaction. The derivation of the cr- 
model from the latter model is straightforward and valid 
for an arbitrary interaction. Unfortunately, the artificial 
supersymmetric model is equivalent to the electron one 
only in the limit of a weak short range interaction. There- 
fore the cr-model we obtained is less general than the 
replica or Keldysh a-models. So, as concerns the pertur- 
bation theory in the diffusion modes or the renormaliza- 
tion group calculations, it is better to perform them with 
the replica or Keldysh cr-models and our goal was not to 
construct a more convenient tool for such calculations. 
Although we have demonstrated how to carry out per- 
turbative and RG computations in the framework of our 
CT-model, the purpose of these calculations was merely to 
check the cr-model. The check is successful and our hope 
is that using the supersymmetric a-model we derived one 
can perform nonperturbative calculations. 

It is well known that for noninteracting systems the su- 
persymmetry technique has real advantages with respect 
to the other methods^ and allows one to consider, e.g., 
localization in wires, level statistics in quantum dots, etc. 
So, our hope is that the supersymmetric cr-model with in- 
teraction derived in the present paper can help us to get 
results in this direction. The supermatrices Q derived 
here should have the correct symmetry (a mixture be- 
tween compact and noncompact sectors) and there is no 
reason to think that, when properly treated, the model 
would lead to wrong results. 

Nevertheless, using the present cr-model is definitely 
more difficult than the conventional cr-model for nonin- 
teracting particles. The problem is that the supermatrix 
Q is now not a8x8or4x4 supermatrix as it was for the 
noninteracting systems. In contrast, the size of the su- 
permatrix Q is now 4M x 4M (for the unitary ensemble), 
where M is the number of the Matsubara frequencies in- 
volved in the calculation (actually, it is infinitely large). 
Therefore, we should be able to calculate integrals over 
the supermatrices of an arbitrary large size. Of course, 
the same problem exists for the replica and Keldysh cr- 
model. In the replica a-model one should, in addition, 
calculate for an arbitrary number N of the replicas and 
set at the end iV = 0. In the Keldysh approach the ma- 
trix Q depends on two energies and one should discretize 
the energy to get a reasonable result. 

The advantage of the supersymmetry approach is that 
the supersymmetric a-model is still simpler and well de- 
fined. 

In principle, the task of making explicit calculations 
with AM X 4M supermatrices does not seem hopeless. 
For calculations in the OD situation one can try to use 
the Itzykson-Zuber integral^SiSi that has been general- 
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ized for the supersymmetric case by Guhr et a\m^. An 
integral with the zero-dimensional supermatrix free en- 
ergy functional taken in the form of Eq. 1)4. 7|) has the 
proper form and one may hope to proceed in this way. 
Another problem is localization in wires. With the su- 
permatrix cr-model one can write using the transfer ma- 
trix technique an effective "Schrodinger equation" and 
then try to solve it. These are the most evident exam- 
ples. Probably, one can try other situations, although 
the proper calculations can be difficult. 

Of course, at the moment we cannot guarantee the 
complete success in this direction. At the same time, it 
is difficult to imagine that either replica or Keldysh a- 
models can be more helpful. Although there is a certain 
progress in reproducing some known results for noninter- 
acting systems using these method£^--4iSiiLi&^ the calcula- 
tions are quite difficult and the results are not complete. 
We do not see how these schemes can be extended to 
interacting systems and believe that the supermatrix cr- 
model developed here is a better opportunity to attack 
the nonperturbative problems. 
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APPENDIX: PARAMETRIZATION OF (/>!, <j>2 
AND (fi, (f2 



Here the definition of the fields (t>i,ipi is given. In 
the following we write supermatrices with 2x2 block 
structure describing spin space, a*-' , b'^ are commuting 
while cr*-' are anticommuting. Variables a, 5, a for differ- 
ent are not related to one another. a(r, r), 6(r, r) 
and a"(r, r), 5"(r, r) are real. We assume to be short 
range and usually set ri(r — r') = ri6{r — r'). Af are 
normalization constants for the Gaussian integrals. Inte- 
gration is over all independent components of the corre- 
sponding matrices. 



and 



(pi(r,n) 



V2(r,«) = 



(/)i(r,n) 



02(r,n) 



a(r, n) 

a(r,n) 



V 







a(r,7i) / 



a}^(y:,Ti) a^^(r,n) 
2^^*(r, — n) a^^(r,n) 





3-^^*(r, — n) a^^(r, n) / 



/ a(r,n) 

a(r,n) 
iT(r, n) 

V a(r,7i) 



cr*(r,-n) 

cr*(r,-n) 

i5(r, n) 

i5(r, n) 



( 

V 



^ — TV" / Dcpi ( ) e~^ / drdr' Str[<pi{r,-n)ur:[^(r-r')<pi(r' ,n)] 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



(A.5) 
(A.6) 
(A.7) 



(•••), 



M j D(f)2 (...) e"* i" Str[02(r,-n),/r-i(r-r')02(r',n)]_ 
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